Magnetoelastic Effects in Iron Telluride 



I. Paul 1 , A. Cano 2 and K. Sengupta 3 
1 Institut Ned, CNRS/UJF, 25 avenue des Martyrs, BP 166, 38042 Grenoble, France. 
2 European Synchrotron Radiation Facility, 6 rue Jules Horowitz, BP 220, 38043 Grenoble, France. 
3 Theoretical Physics Department, Indian Association for the Cultivation of Science, Kolkata-700032, India. 

(Dated: March 9, 2011) 

Iron telluride doped lightly with selenium is known to undergo a first order magneto-structural 
transition before turning superconducting at higher doping. We study the effects of magneto- 
elastic couplings on this transition using symmetry considerations. We find that the magnetic order 
parameters are coupled to the uniform monoclinic strain of the unit cell with one iron per cell, as well 
as to the phonons at high symmetry points of the Brillouin zone. In the magnetic phase the former 
gives rise to monoclinic distortion while the latter induces dimerization of the ferromagnetic iron 
chains due to alternate lengthening and shortening of the nearest-neighbour iron-iron bonds. We 
compare this system with the iron arsenides and propose a microscopic magneto-elastic Hamiltonian 
which is relevant for all the iron based superconductors. We argue that this describes electron-lattice 
coupling in a system where electron-electron interaction is crucial. 

PACS numbers: 74.70.Xa, 74.90. +n, 75.80.+q 



I. INTRODUCTION 

Selenium-doped iron telluride, FeTei-^Se^, is one of 
the several classes of iron based materials which are cur- 
rently being studied for their intriguing superconduct- 
ing properties. 1-5 This material, the so-called 11 system, 
exhibits a magneto-structural transition at low carrier 
doping that is suppressed by the emergence of super- 
conductivity at large enough doping (or under applied 
pressure). 6 This feature of the phase diagram is shared 
by the rest of iron-based superconductors, 7 namely the 
1111 systems with chemical composition i?eOFeAs where 
Re refers to a rare earth metal (e.g., La, Ce and Sm), 8 the 
111 systems with composition ^4FeAs where A refers to 
an alkali metal (e.g., Na), 9 and the 122 systems composed 
of AeFe2As2 where Ae is an alkali earth metal (e.g., Ba, 
Sr and Ca). 10 LiFeAs is the only exception to this rule, 
as the parent compound is already a superconductor. 11 
In Ref. 12 a general phase diagram for the 1111, 122 
and 111 systems was proposed on the basis of symme- 
try considerations which explained the salient features 
of the magnetic and structural transitions. The key in- 
gredient for this general phase diagram is the magneto- 
elastic coupling between the magnetic and the structural 
order parameters. In this paper we show that similar 
symmetry- allowed couplings are also present in the 11 
material, and are important in the sense that they give 
rise to experimentally observable structural distortions 
in the magnetic phase. Our description is based on a 
Ginzburg-Landau mean field analysis of the magneto- 
structural transition of iron telluride followed by a com- 
parison with the iron arsenide systems. Furthermore, we 
propose a microscopic Hamiltonian that provides a uni- 
fied description of the magneto-elastic properties of all 
the iron-based superconductors, and we discuss its im- 
plications. 

Selenium doped iron telluride is represented as 
Fei+ y Tei- x Se x due to the presence of excess Fe in the 



interstitial positions of the Te layer. 13 The amount of ex- 
cess Fe decreases with increasing selenium doping, and 
the phase diagram of the undoped compound is known 
to depend upon the amount of excess Fe. 4 In the follow- 
ing we focus on the system with y w 0.076 in the un- 
doped state (x = 0). It is known to undergo a first-order 
transition at around 65 K from a high temperature para- 
magnetic tetragonal phase to a low temperature antifcr- 
romagnetic monoclinic phase (throughout the paper we 
follow the notations of an unit cell with lFe/cell). 14 The 
transition temperature decreases with increasing x, but 
the structural distortion remains concomitant with the 
magnetic transition. This makes the 11 material differ- 
ent from the 1111 and 111 systems, where there are two 
separate transitions at all doping, and also from the 122 
system where the single transition splits into two tran- 
sitions with doping. More importantly, the symmetry 
breaking associated with the magneto-structural transi- 
tion of the telluride system is different from that of the 
arsenides. Thus, the low-temperature phase of the 11 sys- 
tem (Fig. 1) is monoclinic and is a bicollinear antiferro- 
magnet with wavevector (tt/2,it/2) in the ab plane, 14-16 
in contrast to the low-temperature phase of the iron ar- 
senides (Fig. 2) which is orthorhombic and is a collincar 
antiferromagnet with wavevector (tt, 0). However, de- 
spite these apparent differences, in the following we show 
that a feature common to all the iron based supercon- 
ductors is the presence of important magneto-elastic cou- 
plings whose microscopic origins appear to be the same. 

The magneto-structural transition of Fei+yTei-^Sea, 
has been studied earlier from a field theoretic point of 
view without taking into account all possible magneto- 
elastic terms. 1 At a more microscopic level the 
(7r/2, 7t/2) magnetic structure is somewhat puzzling from 
a weakly correlated band picture because, unlike the case 
of the iron arsenide materials, the magnetic ordering does 
not correspond to a nesting wavevector of the Fermi sur- 
face sheets. This has motivated theorists to seek rationale 
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FIG. 1: Low temperature magneto-structural pattern on the ab plane formed by the Fe atoms of Fei+ H Tei_ x Se x . The spins 
order with wavevector (tt/2,tt, 2) forming a bicollinear antiferromagnet. The associated structural change (c) is a combination 
of a monoclinic distortion (a) of the square lattice, and a dimerization (b) of the ferromagnetic Fe chains due to alternate 
lengthening and shortening of the nearest-neighbour Fe-Fe bonds. 



in models that are based on strong correlation physics 
such as a localized spin model 18 or that describing mag- 
netism induced by orbital ordering. 19 Alternately, it has 
also been argued that doping due to excess Fe changes the 
Fermi surface sufficiently such that (ir/2,ir/2) is indeed 
a nesting wavevector. 20 

The rest of the paper is organized as follows. In the 
next section we describe the magneto-structural transi- 
tion by a Ginzburg-Landau free energy whose form is 
independent of the microscopic details. The (tt/2,it/2) 
magnetic structure is described by four 0(3) magnetic 
order parameters which gives rise to several symmetry al- 
lowed magneto-elastic couplings. Thus, the magnetic or- 
der parameters are coupled to the monoclinic component 
of the uniform strain, as well as to the lattice distortions 
associated with wavevectors of high symmetry in the Bril- 
louin zone. In the magnetic phase the former coupling 
distorts the lattice monoclinically with long and short 
next-nearest-neighbour Fe-Fe bonds (Fig. la), while the 
latter produces long and short nearest-neighbour Fe-Fe 




FIG. 2: Low temperature magneto-structural pattern on the 
ab plane formed by the Fe atoms of the iron arsenide systems. 
The spins order with wavevector (tv, 0) forming a collinear 
antiferromagnet. The associated structural change is an or- 
thorhombic distortion of the square lattice. 



bonds which gives rise to dimerization of the ferromag- 
netic chains (Fig. lb). The experimental observation of 
both these distortions indicate that the magneto-elastic 
terms are non-negligible. 14-16 Furthermore, we find that 
the effective magnetic free energy generated by the lat- 
tice favours a first order transition, as observed exper- 
imentally. We finish section II by briefly recalling the 
mean field description of the iron arsenides, and by com- 
paring it with the iron telluride system. In section III 
we propose a microscopic hamiltonian which compactly 
describes the magneto-elastic properties of all the iron 
based superconductors, and we conclude in section IV by 
pointing out directions for future research. 



II. GINZBURG-LANDAU THEORY 

A. Iron telluride 

The Ginzburg-Landau free energy describing the phase 
transition can be expressed as 

Fgl = Fm + Fe + Fme (1) 

corresponding to the magnetic, elastic and the magneto- 
elastic parts respectively. Hereafter the notation refers 
to the ab plane of the system where the symmetry con- 
siderations are non-trivial. 

The magnetic sector is described by four 0(3) or- 
der parameters which are the Fourier components 
with i = 1, . . . , 4, of the magnetization (M) correspond- 
ing to the wavevectors qi = (ir/2,w/2) = — q 3 , and 
q2 = (— 7r/2,7r/2) = — q4. Thus, at the mean field level 
the magnetization is given by M(r n ) = 53 i=1 exp(«q i - 
r„), where r„ denotes a lattice position. Since M(r„) is 
real, we have Mi = and M 2 = M4. Alternatively, 
one can define four real- valued 0(3) order parameters 
with i = 1, . . . , 4, such that M x = [(Lj + L 3 ) - i(Li - 
L 3 )]/2 and M 2 = [(L 2 + L 4 ) - i(L 2 - L 4 )]/2. In terms 
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of Li the magnetization is 

4 

M ( r «) = X] Li [^(^ ' r «) + sin (q» • r «)l • 

8=1 

The lattice sector is described by means of the strain 
tensor which can be written as €ij(r n ) = Uij + 
I Eq^ofe"; ^) + 9i M i( c l)] eJq r "- Here describes uni- 
form strains and u(q) is the Fourier transform of the 
displacement held u(r„). 21 At the mean held level and 
to the lowest order in an expansion of the free energy 
in terms of the order parameters, we find that the elas- 
tic variables that couple with the magnetic ones are the 
shear component u xy and the quantities Ui = u(qj) with 
i = 5, 6 and 7 corresponding to qs = (7r,7r), q6 = (n, 0) 
and q7 = (0, ir). 

The transformation of the above variables under the 
symmetry operations of the square lattice are summa- 
rized in Table I, using which one can construct Fql from 
the terms that are allowed by symmetry. In the following 
we assume 0(3) symmetry in the magnetic sector, which 
amounts to neglecting spin-orbit coupling. The magnetic 
part of the free energy can be expressed as 



the lowest order magneto-elastic part is given by 

Fme = g\u xy (l\ + l\ — l\ — l 4 ) 

+ g 2 [uf (If -L\-L\ + L\) + u\ [L\ -L\+L\ 
+ g 3 [u% (Li • L 4 L 2 • L 3 ) + u v 7 (L x • L 2 - L 3 • L 4 ; 
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In the above we ignore the standard magneto-striction 
term since it is present in all materials, and is not peculiar 
to the 11 system. 

We minimize Fql with respect to the elastic degrees 
of freedom and we get 
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(Li • L4 — L2 • L3) , 
(Li • L2 — L3 • L4) , 



(5a) 
(5b) 
(5c) 
(5d) 



while the remaining elastic variables are zero at equilib- 
rium. Thus, the lattice mediated effective magnetic free 
energy is given by 



Fm — 



A 



(2) 



where A = A'(T — Tjv) is the only temperature (T) 
dependent coefficient, and Tjy is temperature for the 
magneto-structural transition. The ellipses denote the 
fourth order and the sixth order terms, the latter be- 
ing necessary because eventually the transition is first 
order. At the fourth order there are seven invariant 
terms namely, (i) if, (ii) L\L\ + cyclic terms, (hi) 
(Li • L 2 ) 2 + cyclic terms, (iv) (Li • L 3 )(L 2 • L 4 ), (v) 
(Li ■ L 2 )(L 3 • L 4 ) + 2 o 4, (vi) (Li ■ L 3 ) 2 + (L a • L 4 ) 2 , 
and (vii) L\L\ + L\L\. Thus, at this order there are 
seven independent coupling constants in terms of which 
the phase diagram described by Fm is rather rich. It 
is not the purpose of this paper to investigate how the 
phase diagram varies with different couplings, which is 
anyway a formidable task. The experimentally observed 
magnetic state is described by a non-zero value of any 
one of the Lj, and the energetics of this choice is clearly 
beyond symmetry based arguments and mean field the- 
ory. Consequently, we do not write explicitly the terms 
beyond quadratic order. Next, the elastic part of the free 
energy is given by 
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where cqq is the elastic constant for monoclinic distortion, 
and Sli and Vl 2 are the elastic stiffness of the displace- 
ments at the respective wavevectors. Since the elastic 
sector is not critical, it is sufficient to truncate the ex- 
pansion at the quadratic order. Finally, using Table I, 
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The simple exercise above allows us to make the following 
two points about the experimentally observed magnetic 
phase where L± = L (say) is the spontaneous magne- 
tization and the remaining Li are zero, (i) The lattice 
undergoes a monoclinic distortion with u xy — ~giL 2 /cQQ 
producing long and short next-nearest-neighbour Fe-Fe 
bonds (Fig. la). Simultaneously, the nearest-neighbour 



Fe-Fe bonds dimerize such that 



,y — 



-g 2 L'/n 1 



(Fig. lb). The experimental observation of both these 
distortions imply that the magneto-elastic coupling is 
non-negligible. 14,16 (ii) The lattice mediated effective 
magnetic free energy is given by F' M = —{g\/c§Q + 
<7|/r2i)L 4 /2, which is a negative contribution to the free 
energy at fourth order. As such, this contribution favours 
the experimentally observed first order transition. 



B. Iron arsenides 

To finish this section we briefly recall the magneto- 
elastic properties of the iron arsenides (Fig. 2), studied 
in detail in Ref. 12, in order to compare them with the 11 
material and to facilitate the discussion in the following 



4 



reflections translations 
90° rotation a;-axis x — y axis along x along y 



Li L2 L4 Li L3 L3 

L2 L3 L3 L4 — L4 L4 

L3 L4 L2 L3 — Li — Li 

L4 Li Li L2 L2 — L2 

Uxy Uxy Uxy Uxy Uxy Uxy 

(uf.ug) K,-uf) («S,-ug) («g,«f) (-«S,-«g) (-«S,-«g) 

(uf.tij) («?,-«?) K,-<) («?,«?) (-«!,-«?) OM) 

(uf, u ?) K,-ng) K,-^) Wj) K,^) 



TABLE I: Transformation properties of the magnetic and the elastic order parameters under symmetry operations of the square 
lattice. The definitions of the various order parameters are given following Eq. (1). 



section. The relevant magnetic order parameters for the 
iron arsenides are the Fourier components Mg and M7 
of the magnetization corresponding to the wavevectors 
q6 and q7 respectively. They couple to the orthorhombic 
component u xx — u yy of the uniform strain tensor, and 
at the mean field level the magneto-elastic part has the 
form 22 

F^e = 94 (u xx - u yy ) {Ml - A/ 2 ). (7) 

This term determine (a) if there is a single magneto- 
structural transition (as in the undoped 122 systems) or 
two separate transitions (as in the 1111 or the sufficiently 
doped 122 systems), (b) the order (first versus second) of 
the transitions, and (c) why the system prefers a collincar 
magnetic state instead of a non-collinear order. 

In the case of the iron arsenides, since the system is 
near an orthorhombic transition (this is explicit in the 
phase diagram of the 1111 and the doped 122 systems) 
where the corresponding elastic constant c„ vanishes, it is 
possible to argue that the energy scale g\/\c \ generated 
magneto-elastically is dominant and controls the main 
features of the phase diagram for the magneto-structural 
transition. In contrast, the phase diagram of the 11 sys- 
tem does not exhibit a monoclinic structural transition 
where Cg6 vanishes, and since the magnitude of cqq from 
ultrasound experiments is currently unavailable to us, a 
similar argument for the scale / 1 c 66 1 cannot be made 
at present. Nevertheless, the distortions of the lattice of 
the undoped 122 systems and the 11 systems on enter- 
ing the magnetic phase provide clear evidence that the 
magneto-elastic terms are non- negligible and ubiquitous 
for all the iron based superconductors. 

III. MICROSCOPIC COUPLING 

In the previous section we argued in favour of the ex- 
istence of magneto-elastic couplings in all the iron based 
superconductors. The precise form of the couplings at 
the mean field level differs between the iron arsenides 
and the 11 systems because the magnetic order param- 
eters in these two classes are different. In the following 



we introduce a microscopic magneto-elastic hamiltonian 
whose mean field form captures both Eqs. (4) and (7). 
The purpose of such a hamiltonian is two-fold. Firstly, 
to unify the various magneto-elastic effects observed in 
different classes of systems and provide a common de- 
scription. This is a first step to understand the micro- 
scopic origin of these couplings. Secondly, to go beyond 
mean field theory and study the effect of these couplings 
at the level of fluctuations. 

The simplest microscopic magneto-elastic hamiltonian 
is given by 

HmE = ^2 A<5 ( S " +<5 ' S ") Tn ' S ' \. U ( r n+s) - u ( r ")] ) ( 8 ) 

n , 6 

where n denotes lattice sites, 8 implies nearest- 
neighbours and next-nearest-neighbours, r nt $ — r n+ s — 
r„, \$ is the strength of the coupling that can depend 
on the bond lengths and the bond angles, and S n is the 
electron spin at site n. It is easy to verify that when the 
spins S n and the displacements u(r n ) are replaced by the 
appropriate mean field variables, one obtains the mean 
field results of Eqs. (4) and (7). Thus, while the coupling 
gi = — 4A nnn originates from the next-nearest-neighbour 
magneto-elastic interaction A nnn ; gi-, g-i and 174 are ob- 
tained from nearest-neighbour interaction X nn , implying 
that in the 11 system, in fact, 172 = .93 = — 4A nn . Experi- 
mentally, all the different kinds of distortions in the 1111, 
the 122 and the 11 systems that have been reported are 
such that the ferromagnetic bonds are shorter than the 
corresponding antiferromagnetic bonds, 14,16 ' 23,24 which 
implies that Xs > 0. 

At present it is being debated whether the iron based 
superconductors are weakly interacting metals or those 
that are strongly interacting. In this context it should 
be noted that Hamiltonians of the above type are quite 
common in the study of lattice effects in insulating mag- 
netic systems, where A^ oc dJ/dr is the variation of the 
Heisenberg exchange coupling J with distance. Thus, the 
existence of the magneto-elastic couplings would seem to 
favour the point of view that the magnetic properties of 
these systems are more suitably described by localized 
spin models rather than a weak coupling band picture. 25 
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On the other hand, it is also possible to interpret the 
above Hamiltonian from a band picture with moderately 
strong interaction. Within a semi-classical paramagnon 
theory of magnetism, which is suitable for band met- 
als, S„ can be taken as a paramagnon field which is an 
0(3) variable describing local fluctuations of magneti- 
zation. Therefore, while it is clear that this is an ef- 
fect of electron-lattice coupling in an environment where 
electron-electron interaction is crucial, it is difficult at 
this point to conclude whether a localized limit is neces- 
sary for describing such interaction. 

IV. SUMMARY 

To summarize, in this paper we studied the effects of 
magneto-elastic couplings on the magnetic phase tran- 
sition in Fei+yTei-^Seic. From symmetry considera- 
tions we showed that on entering the magnetic phase, 
these couplings give rise to uniform monoclinic distortion 
(Fig. la) of the unit cell with one Fe/cell, and also induce 
dimerization (Fig. lb) of the ferromagnetic Fe chains due 



to alternate lengthening and shortening of the nearest- 
neighbour Fe-Fe bonds. We also showed that the effec- 
tive magnetic energy generated by the lattice favours a 
first order transition. Finally, we compared this system 
with the iron arsenide systems, and we proposed a mi- 
croscopic Hamiltonian to describe the magneto-elastic ef- 
fects in both these classes of iron based superconductors. 
Microscopically, this is a manifestation of electron- lattice 
coupling in a system where electron-electron interaction 
is strong. In the future we hope to study the effect of 
these couplings on the magnetism and the superconduc- 
tivity. 
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